Journal of Nonlinear Mathematical Physics 



Volume 10, Number 3 f2003). I252W2551 



Letter 



On the Inverse Scattering Approach to the 
Camassa-Holm Equation 

Adrian CONSTANTIN t and Jonatan LENELLS 

Department of Mathematics, Lund University, P.O. Box 118, SE-221 00 Lund, Sweden 
t Corresponding author. E-mail: adrian.constantin@math.lu.se 

Received November 25, 2002; Revised January 06, 2003; Accepted January 07, 2003 

Abstract 

A simple algoritm for the inverse scattering approach to the Camassa-Holm equation 
is presented. 

1 Introduction 

The nonlinear partial differential equation 

m t + 2u x + um x + 2mu x = 0, t > 0, x € K, (1-1) 

in dimensionless space-time variables (x, t) models the unidirectional propagation of two- 
dimensional waves in shallow water over a flat bottom. In (1.1), u(t, x) represents the 
horizontal component of the fluid velocity, or, equivalently, the water's free surface, and 
m = u — u xx is the momentum variable cf . (see also |11| for an alternative derivation) . 
The solitary waves of (1.1) are solitons - they regain their shape and speed after interacting 
nonlinear ly with other solitary waves (see [S] |13j ). While some initial profiles evolve 
into waves of permanent form, others yield waves that break in finite time |3j [I]- Another 
aspect of interest of (1.1) is its bi-Hamiltonian structure ^U] and the induced existence of 
infinitely many conservation laws. This feature is connected to the (formal) integr ability 
of the equation, established in by finding an isospectral problem associated to (1.1). 

Physically relevant cases are solutions of (1.1) with a periodic dependence upon the 
spatial x-variable, as well as solutions which decay at infinity. The objective of this note 
is to prove the integrability of (1.1) for solutions which decay at infinity (the periodic case 
is treated in j^j [7j [5]). For a discussion of the scattering problem for (1.1) we refer to 
[3] and |13j . Recently, in j^j, we proposed an algorithm to solve the inverse scattering 
problem for (1.1). Our aim is to present here a considerable simplification of the approach 
in [H] and to show the applicability of the new approach by means of an example. 
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2 The inverse scattering approach 

If the initial momentum itlq = m(0, •) belongs to the Schwartz class of smooth functions 
/ : R — > R such that for all n%, n 2 > 0, sup \x ni d™ 2 f(x)\ < oo, and m Q + 1 > 0, then both 

properties are preserved by the flow of (1.1) cf. [3] [H]. The isospectral problem in L 2 (R) 
for (1.1) is (see [T]) 

1>xx = -^ + A(m + l)V> (2.1) 

with continuous spectrum (— oo, and at most finitely many eigenvalues in the interval 
(— j,0) cf. P]. The Liouville transformation 



^(y) = ( m (x) + 1 4 ^(x) where y = x + / \/ m (£) + 1—1 

J —oo 



converts (2.1) into 



dy 2 

Here 



+ Q<p = [iip. (2.2) 



4?(y) 4g(y) 16g 2 (y) 4 

with g(y) = m(x) + 1 and spectral parameter \i = — ~ — A. Suitable scattering data for 
(2.1) happen to be the usual scattering data for the problem (2.2), which evolve linearly at 
constant speed under the Camassa-Holm flow cf. jOJ. Therefore the classical Marchenko 
approach (see 0) is applicable and finding Q(t,y) amounts to solving a linear integral 
equation determined by the scattering data for Qo(y), data available from the knowledge 
of mo(x). The only intricate point of this approach is the recovery of m(t, x) from Q(t, y). 
This requires us to solve, given Q, the nonlinear second order differential equation (2.3) 
for q, and then to perform the coordinate transform y i— > x. Equation (2.3) is a Pinney 
equation but the solution for q, given Q, obtained in Jl] is not convenient for our 
purposes (this approach was used in and leads to unnecessary complications). A quite 
intricate (but nevertheless effective) algorithm for the recovery of m was proposed in [£]. 
Below we present an alternative approach that gives a more direct and less complicated 
solution. For convenience we drop the time-dependence in the formulation. 

Theorem Let f(y) be the Jost function at y = oo for the eigenvalue equation 

Vyy = (Q + \) V, (2-4) 
i.e. f is the unique solution of (2.4) with the asymptotic behavior 

f(y) w e~ y/2 and f'(y) ps - - e' y/2 as y oo. 
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If H : R —* R is i/te bijection given by H(y) = J f 2 (£) ' ^ en 

m(x) + l = e 2a / 4 (fl- 1 (e x )) ) s G M. (2.5) 

Proof. Observe that (2.4) is precisely (2.2) with /x = — \. Hence A = 0, and the Liouville 
transformation maps (2.2) into 

4>xx = ^ V>, 

1 

where <p(y) = (jn(x) + 1^ 4 ifi(x). The above equation has the solution e~ x l 2 so that (2.2) 

has the corresponding solution h(y) = (m(x) + l\*e~ x / 2 . But one can easily check that 
the function h has the asymptotic behavior at y = oo prescribed for /, so / = h, i.e. f(y) = 
g 1 / 4 (y) e~ X//2 . This means that = = e x dx. Therefore il(y) = / e s ds = 

/(f) V«) /-oo 

is clearly invertible and y = H (e x ). From the relation m{x) + 1 = q(y) = f 4 (y)e 2x we 
now obtain (2.5). □ 

Remark The previous result reduces the recovery of m(x) from Q(y) to solving the 
linear integral equation 

-y/2, (JH-y)/2 _ Jy-t;)/2 



f(y) = e-W* + J [eW* - e^^J Q{£) /(£) d£, y G R, 

and computing the inverse of the function = / 2 . □ 

J —oo J vsj 

Example A solitary wave for (1.1) is a solution m(t,x) = <3?(x — ct) with a profile f> 
that decays at infinity. Solitary waves can exist only for speeds c > 2, but no mathematical 
expression in closed form seems to be available for $ (see jS]), despite the fact that the 
corresponding potential Q(y) is given explicitly by 

Q(y) = ri = ye 

2c cosh 2 (y-y ) 



with yo £ R (see ^3]). If we choose for simplicity yo = and c = 8/3, then 

1 2cosh 2 ( j/ /4)-l 
Q(y) + I = 8cosh 2 ( y /4) ' yG 



Note that g(y) = 4 cosh 2 (y/4) is a solution to (2.4). But then y \— > p(y) / is also 

/» (?) 

a solution to (2.4). Since / » „ = 5 , we deduce that 

4 g 2 {0 3cosh 3 (y/4) 

/(g) = fi * * , yeR, 

6cosh(y/4) 

and an exact formula for the profile of the solitary wave emerges. □ 
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